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Abstract This chapter emphasizes new magnetic properties that arise when inver- 
sion symmetry is broken in a superconductor There are two aspects that will be 
covered in detail. The first topic encompasses physics related to superconducting 
magnetoelectric effects that arise from broken inversion symmetry. Broken inver- 
sion symmetry allow for Lifshitz invariants in the free energy which can be viewed 
as a coupling between the magnetic induction and the supercurrent. There are simi- 
larities between these invariants and the better known Dzyaloshinskii-Moyira inter- 
action in magnetic systems. These Lifshitz invariants give rise to anomalous mag- 
netic properties as well as new phases in the presence of magnetic fields. Here, we 
will describe the consequences of these Lifshitz invariants, provide estimates for the 
relative magnitudes of the novel effects, and discuss the important role that crystal 
symmetry plays in understanding this physics. Finally, we provide a discussion of 
the fate of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases in broken inversion 
superconductors. In particular, we show how broken inversion symmetry can have a 
profound effect on the stability, existence, and properties of FFLO phases. 



1 Introduction 

One important way in which non-centrosymmetric superconductors differ from con- 
ventional superconductors is in the response to magnetic fields. In particular, the re- 
moval of inversion symmetry leads to new terms in the free energy that give rise to 
magneto-electric effects. These effects are closely related to the appearance of mag- 
netic field generated helical phase in which the superconducting order develops a 
periodic spatial variation. Here we review this physics beginning with a detailed ex- 
amination of the phenomenological theory followed by an overview of microscopic 
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treatments of these problems which include an overview an of the interplay of the 
helical phase and Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases |[I]|2l. 



2 Phenomenology of Single Component Superconductors 

This section reviews the phenomenology relating Lifshitz invariants in the the free 
energy to magnetoelectric effects, vortex structures, and the helical phase. 



2.1 Ginzburg Landau free energy 

A key new feature of non-centrosymmetric superconductors is the existence of Lif- 
shitz invariants in the Ginzburg Landau (GL) free energy |l3]|4]|5]|6]|2l[8l. These give 
rise to magnetoelecti-ic effects 19j[ISl|5][n][l2l[T3], helical phases l6l[T4ll7l[T5llT6ll. 
and novel magnetic properties ||9] |7] [iTl [18] [19] [121 discussed in this chapter. To ex- 
amine the consequences of these invariants we initially consider a GL theory for a 
single component order parameter (for example, an 5-wave superconductor) and add 
the most general Lifshitz invariant allowed by broken inversion symmetry. Specific 
Lifshitz invariants are tabulated in Table 1 for different point group symmetries of 
the material in question. Since the primary goal is to reveal the new physics arising 
from these invariants, we ignore the role of any anisotropy that might appear in the 
usual GL free energy. Under these conditions the GL free energy under considera- 
tion is (we work in units such that h = c= 1): 

j d\\a\^\^+K^*D^^+ K,jB, [,7 * (d, ,7 ) + 77 (D.-t] )*] + | | 77 | V ^ | , 

(1) 

where a = <X(){T — 7^), D, = —/V, — 2M, and B = V x A. From this free energy, 
the GL equations can be found by varying the above with respect to A and T7 . This 
results in the following: 

ari+li\ri\^ri + KD^ri+Kij[2hi{Djri) + iriVjBi] = (2) 

and 

J/ = ^[Vx (B-47rM)],--2e/:[77*(A-77) + 77(A-i7)1+4e/:,-iB,|T7p (3) 

where 

M, = -/:,-,[77*(D,77) + T7(D,77)*]. (4) 

These equations are joined by the boundary conditions (which follow from the sur- 
face terms that arise from integration by parts in the variation of F): 
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[Kni{Diri)+KijBifijri]houndary = (5) 

where fij is the component of the surface normal along j, and the usual Maxwell 
boundary conditions on the continuity of the normal component of B and the trans- 
verse components of H = B — 4-nM (the appearance of M due to the Lifschitz in- 
variants makes this boundary condition non-trivial). Note that adding the complex 
conjugate of Eq.|5]multiplied by rf* to Eq.|5]multipUed by 77 yields J- «|;,o„„jarj = 0. 



Point Group 


Lifshitz Invariants 





K{BJ, + B,j, + BJ,) 




T 


KiBJ, + B,j, + BJ,) 




Db 


Ki {BJ,- + Byj, + B-J,) + K2BJ, 




Cev 


K{B,j,-B,h) 




C6 


Kt (S,/,- + B,j, + BJ^ + K2BJ, + Ki {BJ, 


-By id 


Da 


Ki [BJ, + B,j, + BJ^ + K2BJ, 




Cav 


KiBJ,-B,j,) 




Did 


K{B_,j,-B,j,) 




Ca 


Ki (B,/,- + B,jy + BJ^ + K2BJ, + Ki {BJ, 


-Byj.y) 


Sa 


Ki [BJ, - Byjy) + K2 (B,;, + B,y,) 




Dt. 


Ki {BJ,- + B,j, + BJ^ + K2BJ, 




Cir 


K{BJ,-B,j,) 




C3 


Ki {BJ,- + BJ, + BJ^ + K2BJ, + {BJ, 


-Byj,) 


O2 


K^BJ, + K2BJ, + KiBJ, 




C2„ 


KiBJ, + K2B,j, 




C2 


KiBJ, + KiByky + Kj,BJ, + KaB,], + K5BJ, 


Cv 


KiB.k, + KiBJj + K3BJ, + KAByj, 




Cl 


all components allowed 





Table 1 Allowed Lifshitz invariants for different point groups. Here y, = r]*(D,r]) + 77(0,7])*. 

The appearance of M in Eq. |4] and the associated magnetization current leads 
to new physics in non-centrosymmetric superconductors. Also note, as is the case 
for centrosymmetric superconductors, the boundary conditions are valid on a length 
scale greater that £,0, the zero-temperature coherence length. In the following few 
subsections, we present the solution to some common problems to provide insight 
into the role of the Lifshitz invariants. 



2.2 Solution with a spatially uniform Magnetic field: Helical Phase 

In situations when the magnetic field is spatially uniform, the GL equations describ- 
ing the physics can be greatly simplified by introducing the following new order 
parameter; 

ri = riexp liq-x) = 77 exp U—^ — ) • (6) 
The GL free energy for fj no longer has any Lifshitz invariants and is 
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a-B,K,,nBjKj„, \f,\^+Kifi*T)-fi + !^\fif + —\. (7) 




The resulting new GL equations are now those of a single component superconduc- 
tor with a magnetic field induced enhancement of (this magnetic field enhance- 
ment is discussed in more detail in Chapter 1). These new GL equations follow from 
a minimization of Eq.|2]with respect to A and fj. Note that the phase factor intro- 
duced above cancels the additional current contribution from the Lifshitz invariants 
in Eq. [5] and also cancels the related Lifshitz invariant contribution to the bound- 
ary condition. Furthermore, the magnetization that follows from Eq.|7]by taking the 
derivative with respect to B, coincides with that due Eq.|4]found prior to the redefini- 
tion of the order parameter This modified free energy of Eq. ^immediately implies 
that some results from the usual GL theory apply. In particular: 

i) the vortex lattice solution near the upper critical field is the same as that of 
Abrikosov. 

ii) the surface critical field //^j is the same as that of DeGennes. The order cor- 
rections to Tc do not change //^3 to leading order in {T^ — T) /Te- 
rn) the critical current in this wires will show no unusual asymmetry (this conclu- 
sion differs from that of Ref. I?]). 



2.2.1 Helical Phase 

The main new feature that appears in a uniform magnetic field is the spatial mod- 
ulation of the order parameter Since 7] develops a helical spatial dependence in 
the complex plane, the resulting thermodynamic phase has been named the helical 
phase. Since helicity of the order parameter is related to its phase, an interference 
experiment based on the Josephson effect would provide the most reliable test to ob- 
serve this. Indeed, such an experiment has been proposed Q . In particular, consider 
the example of a 2D non-centrosymmetric superconductor (with a Rashba spin-orbit 
interaction) with a Zeeman field applied in the 2D plane. Then consider a Josephson 
junction between this and another thin film superconductor that is centrosymmetric. 
For a magnetic field applied in the plane of the film perpendicular to the junction 
and with the non-centrosymmetric superconductor oriented so that the helicity q is 
perpendicular to the field ; we find this gives rise to an interference effect analogous 
to the standard Fraunhofer pattern. For this experiment, the film must be sufficiently 
thin so that the magnetic field and the magnitude of the order parameter are spatially 
uniform. 

To illustrate this, consider the following free energy of the junction 




(8) 



where the integral is along the junction. The resulting Josephson current is 
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Ij = Im 



(9) 



Setting the junction length equal to 2L, and integrating yields a maximum Josephson 
current of 



Ij = 2t\1'l'\\'Fi' 



|sin(gL)| 
\qL\ 



(10) 



This demonstrates that the Josephson current will display an interference pattern 
for a field perpendicular to the junction. Note that in the usual case the Fraunhofer 
pattern would be observed for a magnetic field perpendicular to the thin film for 
which a finite flux passes through the junction. 



2.2.2 Magnetoelectric Effect 

Amongst the early theoretical studies of non-centrosymmetric superconductors, it 
was pointed out that a supercurrent must be accompanied by a spin polarization of 
the carriers ifTOl . Within the macroscopic theory given above, this spin polarization 
is described by the magnetization in Eq. |4] This magnetization appears when the 
supercurrent is non-vanishing due to a finite phase gradient. Subsequent to this pro- 
posal, it was suggested that the converse effect would also appear: a Zeeman field 
would induce a supercurrent [[5]. This would follow from the expression for the cur- 
rent of Eq.|3]when the usual GL current (2eK[ri*{Diri) + T7(D,T7)*]) vanishes. How- 
ever, the latter proposal does not include the possibility discussed above that the 
order parameter develops a spatial modulation in the presence of a spatially homo- 
geneous magnetic field (which leads to a nonvanishing 2eK[ri* (Dir]) + 77(0,77)*]). 
Indeed, this new equilibrium state ensures that the resultant supercurrent is vanish- 
ing. Nevertheless, as pointed out in Ref. llT4l . it is possible to create this current 
using a geometry similar to that used to observe Little-Parks oscillations. In par- 
ticular, the supercurrent has two contributions, one is the current due to the Lifshitz 
invariants and the other is the usual GL current 2e/r[T7*(D,T7) + 77 (D;T7)*]. In the he- 
lical phase, these two contributions exactly cancel. By wrapping the superconductor 
in a cylinder, the condition that the order parameter is single valued does not allow 
the helical phase to fully develop since arbitrary spatial oscillations are not allowed. 
Consequently, when a magnetic field is applied along the cylindrical axis, a non- 
zero current can flow. The resulting current will develop a periodic dependence on 
the applied magnetic field lfT4l . 



2.3 London Theory and Meissner State 

We now turn to situations in which the magnetic field is not spatially uniform. The 
Lifshitz invariants lead to new physics for both the single vortex solution and for the 
usual penetration depth problem. To see this, we begin with the London limit and 
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set 7] — \ri\e and assume that the magnitude |t]| is fixed. The GL free energy is 
then minimized with respect to 9 and A. The minimization with respect to yields 

KiW ■{Ve-2eA)+KijViBj^0 (11) 

which is equivalent to the continuity equation for the current (V • J = 0). The mini- 
mization with respect to A yields 

J' = ^ (B -47rM)], = -^[A, - ^V,0 -^(T,-,-B,] (12) 



J 



with 



4nMi = -L^cj, .(A,- - lv,0), (13) 



2e 

J 



1 






(14) 


1 






(15) 




(16) 



1/A^ = ^n{2e)'^K\ri\^ and (7,y = 16neX^Kij. We take the surface normal is along the 
Z direction and that the applied field is oriented along the y direction. Note that by 
applying an appropriate rotation to the fields in the free energy, this geometry results 
in no loss of generality. We assume that there are spatial variations only along the 
direction of the surface normal (z). We therefore have from V • B = that B- = 0. We 
further choose A = [A_v(z),Av(z),0] sothatB = {—dAy/dz,dAx/dz,0) and work in 
a gauge where VO = 0. The three components of Eq.[T2lvields 

dBy I d ^ ,1 

Note that contributions from a^y and (7,., cancel in the above. Taking derivatives of 
Eq . fT4l and [TS] with respect to z, using Eq.[T6]to eliminate A-, we find 

ald^By_ 1 a,, + a,,dB, a,ya-,d^B^ 
^ l^' dz^ ' A2 dz A2 dz^ ^ ' 

^ X'^' dz^ ^ V^ ' A2 dz A2 dz^ ■ ^ 
The above must be solved with the boundary conditions B, (z °°) =0 and 

Hy= By{z = 0)-4KMy{z = 0) (19) 
0= By{z = 0)-AnM^{z = 0) (20) 

where Hy is the applied field. Mx,My can be found using Eq.[T3]and Eq.[T4l[T5] and 
[161 to ehminate Av,Ay, and A, in favor of and By and their derivatives. By setting 
Bi — B,oexp( — 5z/A), the solution can be found analytically. The general form of 
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2.3.1 O point group 

A representative material is Li2Pt3B 1201 1211 . This problem has been solved in 
Refs. ll9l [T2II . In this case there is only one Lifshitz invariant: KiB • j. Since this 
is a scalar under rotations the solution is the same for any orientation of the surface 
normal. The equations for B become; 

d^By 1 8 <9B, 
d^B, _ I 8_dB^ 

~d^^ A2 ■'^"F^- ^^^^ 

where 5 = — 2(J,v( (note Gxx = Cfyv in this case). This coupled set of equations can 
solved for B±= B^ ± iBy ||9] [l2l with the result that to first order in 8/X: 



5 5 3 

By= //,[cos^ + -sin^]e-'/^ (23) 



//v[|cos||-sin||]e-^V^. (24) 

Physically, this implies that the the magnitude of the is discontinuous as it crosses 
the surface (though not that of Z?,,) and that B also rotates inside the superconductor 
Note that in a slab geometry, is of opposite sign on the two sides of the slab. It 
may be possible to observe this through muon spin resonance experiments. 



2.3.2 C4v point group 

A representative material is CePtaSi ll22l . In this case, the single Lifshitz invariant 
is generated by a Rashba spin-orbit coupling and is given by K\z-B x j. This im- 
plies a = Oxy = — C7yT ^ 0. The solution of the London problem now depends upon 
surface orientation and has been considered in Ref . llTSI . We consider two situations 
here: the surface normal along and perpendicular to z (the four-fold symmetry axis). 
Consider first the normal along the z direction (in this case the applied field is Hy 
and we find that B, = 0), then we have the usual London equation 

d^By 1 



with the unusual boundary condition Hy\-=Q = {By + jBy)\^=o. This yields the so- 
lution 
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Byiz) = TT^^-'^' (26) 

These equations show that there is no rotation of B across the sample surface. How- 
ever, the magnetic induction B is discontinuous as the surface is crossed. Again, in 
a slab geometry, the discontinuity in By is opposite for the two sides of the slab. 

For the surface normal perpendicular to the z direction, the situation is different. 
To be concrete, consider the normal along the i direction and the applied field along 
the y direction (for the field along the z direction the usual London Equations result). 
In this case, it is again permissible to set = and solve for By to find 

By = -^e-'/^ (27) 

~ 2 

where a — Oxy and — p-). 



2.4 Spatial structure of a single vortex 

The London theory can also be used to examine the field distribution of a vortex in 
a strongly type II superconductor Again, the lack of inversion symmetry introduce 
some new physics. Here we focus (as above) on two examples with point groups 
O and C4v and provide the solutions of Refs. ifTsl [191 [121 [T3l . The approach used 
in these publications is to consider the parameter Gij/X to be small and then the 
Lifshitz invariants perturb the usual London solution. When there are no Lifshitz 
invariants, the solution to the London equations are = —0 (0 is the polar angle) 
and the field is appUed along the « direction 

B = ^Koir/Xyi (28) 

where Ko{x) is a modified Bessel function. The perturbative solutions depend upon 
the specific form of the Lifshitz invariants and we turn to a discussion of two cases 
in turn. 



2.4.1 O point group 

The solution in this case was found in Ref. |[T2l[T3ll . The modified London equation 
is (the problem does not depend upon field direction) 

The new term implies that, in addition to the field along z, there is an additional 
component along ^ . The authors of Ref. HJ] [13] find that to first order in d /X the 
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additional field is 

'W^i(^) (30) 
where /i and Ki are modified Bessel functions of the first kind. 



2.4.2 point group 

The solution in this case was found in Ref. lfT8l[T9l . The fields that appear due to 
the Lifschitz invariants depend in this case upon the orientation of the field. For the 
field along the y direction, it is found that the solution for B is given by (correct to 
first order in a/A) lfT8l[T9ll 

^-^2mr+jm)y- (31) 

Physically, this impUes that the maximum value of By is shifted from the vortex 
center This shift has also been seen in a full numerical solution of the Ginzburg 
Landau equations ITtI . For the field along the z direction (the four-fold symmetry 
axis), the B field is unchanged and there is an induced magnetization along the radial 
direction (TSl (this radial magnetization was also found in the vortex lattice solution 
mar Ha Hi)- 



2.5 Vortex Lattice Solutions 

For fields near the upper critical field, there have been a variety of studies on the 
Abrikosov vortex lattice ||7]|24]|23]|25]. Some of these studies predict multiple phase 
transitions in the vortex lattice state ll24l l23l l25l . These studies are based on micro- 
scopic weak-coupling theories and involve an interplay of paramagnetism, orbital 
diamagnetism, gap symmetry, band structure, and spin-orbit coupling 1241 l23l |25]| . 
While this chapter will not address these vortex lattice transitions, we will address 
some of the microscopic issues in the next chapter Here we focuss on the GL theory, 
for which the predictions are more straightforward. In particular, near the upper crit- 
ical field, the magnetic field is approximately uniform and the considerations above 
imply that the vortex lattice is hexagonal (perhaps distorted by uniaxial anisotropy). 
Consequently (following the arguments of Section II B) , the order parameter solu- 
tion near the upper critical field is 77 (r) = c«5f exp(/q • r)<pQ{x,y) where <pQ{x,y) is a 
lowest Landau level (LLL) solution. This solution, combining a phase factor and a 
(LLL) solution, has been called the helical vortex phase. The primary consequence 
of this solution is that the upper critical field is enhanced due to the presence of 
the Lifshitz invariants Q . We note that due to the degeneracy of the LLL solution. 
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there is ambiguity in the existence of the phase factor. In particular, the LLL so- 
lution ^Q{x,y) = e''^y''^'H(j)Q(^x,y — Ty) (/// is the magnetic length) is degenerate with 
^(){,x,y), consequently in some circumstances the wavevector q can be removed in 
favor of a shift of origin. This can be done whenever q is perpendicular to the ap- 
plied magnetic field (this is the case for €4^ point group symmetry but not for O 
point group symmetry). We feel that is still meaningful to speak of the helical vor- 
tex phase for the point group €4^ because the same phase factor implies an increase 
of the in-plane critical field in two-dimensions for which this ambiguity does not 
exist. The name helical vortex phase reveals the link between the solutions in two 
and three dimensions. 

In addition to studies near the upper critical field, there has been one numerical 
study of the time dependent GL equations in the vortex phase fTTl . This study found 
the surprising result that the vortices flow spontaneously, in spite of the lack of an 
applied current. The claim is that the paramagnetic supercurrent (the magnetization 
current V x M) is the origin of this spontaneous flux flow. We note that in this study 
the following boundary condition was used; Routside = ^imide- This differs from the 
continuity of H = B — 4-nM discussed above. In the problem that was studied, M is 
non-trivial and an examination of its neglect in the boundary condition can be seen 
to be equivalent to having a current flow. We argue that this current is cause the 
spontaneous flux flow. We note that the boundary conditions discussed here should 
be used in problems where the minimum length scale is ^0, the zero temperature 
coherence length. However, at lengths scale smaller than this, a microscopic theory 
is required and the single particle quantum mechanical wavefunctions will obey 
quite different boundary conditions. 



2.6 Multi-component order parameters 

There have not been as many studies on Lifshitz invariants in non-centrosymmetric 
superconductors in cases when the order parameter contains more than one complex 
degree of freedom. There has been one noteworthy result, which is the appearance 
of the helical phase when no magnetic fields are applied 1211 [8ll. In particular, if the 
ground state of the multi-component order parameter breaks time-reversal symmetry 
II26II27I . then the lack of both parity and time-reversal symmetries allows the helical 
phase to appear. As an example, consider the three dimensional irreducible repre- 
sentation of the point group O, with an order parameter 77 where the components 
transform as the {x,y,z) component of a vector. The following Lifschitz invariant 
exists ISl 

iK{ri*iDyri3 + JTIA.??! + 173^x^72 - c.c). (32) 

This Lifschitz invariant leads to a ground state order parameter 77 = e"''(l,/,0). The 
state 77 = (1,!,0) breaks time reversal symmetry and thus mimics the role of the 
magnetic field in the single component case. 
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The phenomenological arguments of the previous section have also been the subject 
of many microscopic calculations. These calculations, while all related, focus on and 
extend different aspects of the phenomenological theory above. In particular, four 
points of contact exist between the phenomenological theories and the microscopic 
theories. These are: direct calculations of the Lifshitz invariants in the free energy 
in Eq. [T] calculations of the magnetization in Eq. |4l calculations of the current in 
Eq. [3j and calculations of the helical wavevector q in Eq. |6] We briefly review the 
first three of these and then turn to a more complete overview of microscopic studies 
of the helical phase since this turns out to be closely linked to FFLO phases. 



3.1 Contact between microscopic and macroscopic theories: 
Lifshitz Invariants 

The direct calculation of the Lifschitz invariants in Eq.[T]has been carried out by a 
few authors Q |4l |5] [D and can be found in Chapter 1 of this book. In particular, the 
non-interacting Hamiltonian is 

Ho^Y. L [^(k)5„,3 + 7(k)-(T„^]«L«k/3 (33) 

where a^^, (aka) creates (annihilates) an electronic state |ka), ^(k) = £(k) — jj, 
denotes the spin-independent part of the spectrum measured relative to the chemical 
potential ji, a,l5 =tii ^6 spin indices, a are the Pauli matrices, and the sum over 
k is restricted to the first Brillouin zone. In the helicity basis, this Hamiltonian is 
diagonaUzed with energy bands given by 

^±(k) = ^(k)±|7(k)| (34) 

with the Hamiltonian 

^o = E L ^>L(k)4'^kA, (35) 

k l=± 

where the two sets of electronic operators are connected by a unitary transformation, 

flka =X!«aA(k)ckA, (36) 

with 

(«u(k), Mu(k)) = . i/i I , ^ f — ■ ^^^^ 

In the limit that only one of the bands cross the the Fermi energy (this can be realized 
for superconductivity at the surface of a topological insulator [l28]| ), the following 
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weak-coupling result for the coefficients defining the Lifishitz invariants of Eq.[T]is 
found 

^,; = -^^(0^(k)1^-(k)v,,-(k)) (38) 

where A'o is the density of states of the band at the chemical potential, 0(k) de- 
scribes the superconducting state and is an even function belonging to one of one- 
dimensional representations of the point group of the crystal, (...) means the aver- 
aging over the Fermi surface, jis is the Bohr magneton, and 

«n^-ri^.2|p>. ,39, 

Eq. [38] is valid when there is only a single band present. When two bands are 
present (as is often the case), and assuming that (j){k) is the same for both bands, 
then Eq. [38]must be multiplied by the factor 

5N = {N+-N^)/{N+ +N^). (40) 

where N± are the density of states of the two bands (A^o = +N-). Microscopic 
calculations of the Lifshitz invariants are limited to the regime near Tc where the GL 
theory is valid. 



3.2 Contact between microscopic and macroscopic theories: 
current and magnetization 

In the limit of small magnetic fields (B) and small phase gradients (V0) in the super- 
conducting order parameter, it it possible to find microscopic extensions to Eq.[3]and 
Eq. |4]that are valid for all temperatures. This has been carried out in Refs.llSl fTOll . 
Here, we follow the notation of Ref. ijS). In the clean limit, for 2D cylindrical bands 
with a Rashba interaction (/(k) = an x p(k)) Eq.|3]and Eq.|4]can be rewritten as 

Jx — P.V 2m '^°y 

My = f riV.,0 (41) 

where My is the magnetic moment, is the superfluid density, and 

'^(^) = -r72\l^F+{^-Y{T,A+}~pF^{\-Y{T,A^)}] (42) 

where pf ± are the Fermi momenta for the two bands, A± are the gaps on the two 
bands, /i is the Fermi energy, and Y{x) is the Yoshida function. Note that Eq.l42]is 
proportional to 5N in the limit 5N < < 1 . 

The role of Fermi liquid corrections has also been examined [TTI| in this con- 
text. This study has found the that the only Fermi liquid corrections that alter the 
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current contribution from the Lifshitz invariants are ferromagnetic correlations. If 
there are no ferromagentic correlations, then Eq. l42l is unchanged. This is important 
in heavy Fermion materials, where the effective mass enhancement suppresses the 
usual supercurrent but does not change Eq. l42]|iin . 



3.3 Microscopic Theory of the Helical and FFLO Phases 

The helical phase has received a great deal of attention from the microscopic point 
of view Il29l|7l|30l|25l|23|25l[ll[l5l[l6l|28l. One reason for this is that it is 
closely related to the FFLO phase HI |2l in which the superconducting order param- 
eter develops a periodic spatial structure. The interplay between these two phases 
is not trivial. It is perhaps not surprising that spatially oscillating superconductor 
solutions readily appear in non-centrosymmetric superconductors when magnetic 
fields are applied. In particular, a state with momentum k at the Fermi surface will 
generally not have a degenerate partner at — k with which to form a Cooper pair 
when both parity and time reversal symmetries are broken. The state k would rather 
pair with a degenerate state — k + q and in this way generate a spatially oscillating 
superconducting order parameter 

Fig. 1 A magnetic field di- 
rected as shown together with 
a Rashba spin-orbit interac- 
tion shifts the center of the 
large and small Fermi surfaces 
by ±q/2. The sinaller dot rep- 
resent the point (0, 0) (center 
of Fermi surfaces without 
field) and the two larger dots 
represent the points (0, —q/2) 
and (0, q/2) (centers of the 
new Fermi surfaces). Pair- 
ing occurs between states 
of k + q/2 and -k + q/2, 
leading to a gap function 
that has a spatial variation 
A[x) = /loexp(iq-x). From 

Ref. nn. 



The microscopic origin of the spatially oscillating states can be understood by 
an examination of the single particle eigenstates when a Zeeman field H is included 
(for now we ignore the vector potential A) 




k,a,/3 



(43) 
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The single particle excitations now become 

^± (k, H) = ^ (k) ± y^72(k)-2AiB7(k)-H + AijH2. (44) 

In the limit |7| >> \H\, this becomes (we ignore the small regions of phase space 
for which 7=0) 

^±(k,H)«^(k)±MB7(k)-H. (45) 

The origin of pairing states with non-zero q (that is A (x) e'l") follow from this 
expression. As an example, consider a Rashba interaction 7 = 7x {kyX — k^y) for 
a cylindrical Fermi surface and a magnetic field along x. In this case, as shown 
in Fig. [T] the Fermi surfaces remain circular and the centers are shifted along the 
y direction. A finite center of mass momentum Cooper pair is stable because the 
same momentum vector q can be used to pair every state on one of the two Fermi 
surfaces. In the more general case, for a non-zero q state to be stable, the paired 
states should be degenerate: ^±(k + q,H) = k + q,H), this gives the condition 
hq-\f = ■ 7(k). This differs from the condition for the usual FFLO phase, 
for which hq-\f = /iB|H|. The optimal paring state corresponds to finding q that 
satisfies the pairing condition for the largest possible region on the Fermi surface. 

The above paragraph also reveals the origin of the interplay between the helical 
and FFLO phases. In particular, the two Fermi surface sheets prefer pairing states 
with opposite sign of q. Choosing a particular q allows pairing on one Fermi surface, 
but not on the other. This naturally leads to competition between single-q (helical) 
and multiple-q (FFLO-like) states. Which state appears depends upon the details 
of the system. Without going into further microscopic details, which can be found 
in Refs. ||29l|21[30l|25l|23l[ll[l5l[l6l |2l, we summarize some of the main re- 
sults here. One important result is that since there are two sources of the modulation 
q (FFLO-like physics and Lifshitz invariants), there are two typical values for the 
magnitude of q |fT5l[T6ll30ll251 that both appear in different regions of the tempera- 
ture/magnetic field phase diagram. In particular q w H^b/vf stems from FFLO-like 
physics related to Fig.[T]and is the value of q in the high-field regime (in clean ma- 
terials). While q « 5NHij.b/vf stems from the Lifschitz invariants and is the typical 
magnitude of q in the low-field regime ifTSl [161 [30l . As shown in Fig.|2] in the clean 
limit, both single-q and multiple-q phases exist lfT6l [TSl . However, the multiple-q 
phase become less stable as 5N increases lfT6l . We note that in the case of super- 
conductivity at the surface of a topological insulator, which is akin to 5N = 1, only 
the single-q exists (281 . In the dirty limit the multiple-q phases no longer appear, 
while the single-q phase with q w SNHfis/vF is robust ll30l ITSl . Finally we note 
that when the vector potential is also included then novel vortices and vortex phases 
may appear IB] [3l] llll |24l |25l . 



Magnetoelectric effects, helical phases, and FFLO phases 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



T/Tc T/Tc 

Fig. 2 Typical phase diagram showing both multiple-q and single-q (helical phase) phases as a 
function of Zeeman field in a clean non-centrosymmetric superconductor for two different values 
of 5N. These calculations where canied out with a Rashba spin-orbit interaction and a 3D spherical 
Fermi surface (a 2D cylindrical Fermi surface gives similar results). For fields Hlig/Tc < 1.5, 
q Ki SNHjlg/vf (for the SN = 0, this leads to q = 0), while for higher fields q ss Hjlg/vf. From 
Ref. [161 . 



4 Conclusions 

In this chapter we have examined the role of Lifshitz invariants that appear in the 
Ginzburg Landau free energy of non-centrosymmetric superconductors. These in- 
variants lead to magnetoelectric effects, novel London physics in the Meissner state, 
new structure in individual vortices, and a helical phase in which the order parameter 
develops a periodic spatial variation. Additionally, we have provided an overview of 
theoretical developments in the microscopic description of this physics. 
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